Superfluid behaviour of a two-dimensional Bose gas 
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Two-dimensional (2D) systems play a special role in many-body physics. Because of thermal 
fluctuations, they cannot undergo a conventional phase transition associated to the breaking of a 
continuous symmetry pQ. Nevertheless they may exhibit a phase transition to a state with quasi- long 
range order via the Berezinskii-Kosterlitz-Thouless (BKT) mechanism A paradigm example is 
the 2D Bose fluid, such as a liquid helium film 3 , which cannot Bose-condense at non-zero temper- 
ature although it becomes superfluid above a critical phase space density. Ultracold atomic gases 
constitute versatile systems in which the 2D quasi-long range coherence and the microscopic nature 
of the BKT transition were recently explored 4- 6 . However, a direct observation of superfluidity in 
terms of frictionless flow is still missing for these systems. Here we probe the superfluidity of a 2D 
trapped Bose gas with a moving obstacle formed by a micron-sized laser beam. We find a dramatic 
variation of the response of the fluid, depending on its degree of degeneracy at the obstacle location. 
In particular we do not observe any significant heating in the central, highly degenerate region if 
the velocity of the obstacle is below a critical value. 



'Flow without friction' is a hallmark of superfluidity 
[7]. It corresponds to a metastable state in which the 
fluid has a non-zero relative velocity v with respect to 
an external body such as the wall of the container or an 
impurity. This metastable state is separated from the 
equilibrium state of the system (v — 0) by a large energy 
barrier, so that the flow can persist for a macroscopic 
time. The height of the barrier decreases as v increases, 
and eventually passes below a threshold (proportional to 
the thermal energy) for a critical velocity v c . The micro- 
scopic mechanism limiting the barrier height depends on 
the nature of the defect and is associated to the creation 
of phonons and/or vortices [7]. While the quantitative 
comparison between experiments and theory is compli- 
cated for liquid 4 He, cold atomic gases in the weakly in- 
teracting regime are well suited for precise tests of many- 
body physics. In particular, superfluidity was observed 
in 3D atomic gases by stirring a laser beam or an optical 
lattice through bosonic or fermionic [12] fluids and 

by observing the resulting heating or excitations. Here 
we transpose this search for dissipation-less motion to a 
disc-shaped, non homogeneous 2D Bose gas. We use a 
small obstacle to locally perturb the system. The obsta- 
cle moves at constant velocity on a circle centred on the 
cloud, allowing us to probe the gas at a fixed density. We 
repeat the experiment for various atom numbers, temper- 
atures and stirring radii and identify a critical point for 
superfluid behaviour. 

Our experiments are performed with 2D Bose gases 
of N = 35000 to 95000 87 Rb atoms confined in a cylin- 
drically symmetric harmonic potential V(r) + W(z) [14] . 
The trap frequencies are w r /27r = 25.0(5) Hz in the hor- 
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FIG. 1: Stirring a 2D Bose gas. a, A trapped 2D gas 
of 87 Rb atoms is perturbed by a focussed laser beam, which 
moves at constant velocity on a circle centred on the cloud. 
The stirring beam has a frequency larger than the 87 Rb reso- 
nance frequency ('blue detuning' of w 2nm) and thus creates 
a repulsive potential which causes a dip in the density pro- 
file, b, The stirring beam is focussed onto the 2D cloud via 
a microscope objective of numerical aperture 0.45, which is 
also used for imaging. We overlap the two beam paths with a 
polarizing beam splitter cube (PBS). The position of the stir- 
ring beam is controlled by a two-axis piezo-driven mirror, c, 
in situ false-color image of the 2D cloud in the presence of the 
laser beam (average over six images). From the dip in the den- 
sity we deduce the waist of the laser beam as Wo = 2.0(5) /im. 
In this image, the intensity of the beam is chosen three times 
higher than in the stirring experiment to make the hole well 
visible even in the center of the cloud. We use similar im- 
ages, but with the stirring beam switched off, to determine 
the temperature T and the chemical potential fj, from a fit of 
the Hartree-Fock prediction to the wings of the cloud [14] , 
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FIG. 2: Evidence for a critical velocity. Two typical curves of the temperature after stirring the laser beam at varying 
velocities, a, In the superfluid regime, we observe a critical velocity (here v c = 0.87(9) mm/s), below which there is no 
dissipation, b, In the normal regime, the heating is quadratic in the velocity. The experimental parameters are (TV, T, fi, 
r) = (87000, 89 nK, k B x 59 nK, 14.4 //m) and (38000, 67 nK, k B X 39 nK, 16.6 /im) for a and b, respectively, yielding fi l0 c/k B T = 
0.36 and jJ.\ c/k B T = 0.04. The data points are the average of typically ten shots. The y error bars show the standard deviation. 
The x error bar denotes the spread of velocities along the size of the stirring beam (1/^/e radius). The solid line is a fit to 
the data according to equation |T]). Note that the three low-lying data points in a correspond to the completion of half a turn 
and are correlated to a displacement of the cloud, which may be responsible for the observed 1.5 nK temperature shift of these 
points, c and d, Calculated radial density distribution for the clouds in a and b, respectively. The dashed blue curve shows the 
superfluid density, the solid red curve shows the normal density. The stirring beam potential is indicated by the grey shaded 
area (in arbitrary units). The densities are calculated via the local density approximation from the prediction for an infinite 
uniform system [15]. The jump of the superfluid density from zero to a universal value of 4/Aj? 1B (where AdB is the thermal de 
Broglie wavelength) is a prominent feature of the BKT transition. The normal density makes a corresponding jump to keep 
the total density continuous. 



izontal plane and lo z /2tt = 1.4 (1) kHz in the vertical di- 
rection. We use gases with temperature T and central 
chemical potential \i in the range 65-120 nK and /cb X (35- 
60) nK, respectively. The interaction energy per particle 
is given by [/ int — (h 2 g/m)n [I3] . where n is the 2D 
spatial density (typically 100 atoms//zm 2 in the center), 
m the atomic mass and g the dimensionless interaction 
strength. Here g = ^/&na/l z = 0.093, where a = 5.3 nm 
is the 3D scattering length and l z = yjh/mu) z [Mj. The 
energy hui z (fcs x 70 nK) is comparable to k^,T and U- m t 
(~ k-g, x 40 nK at the trap center), and the gas is in the 
quasi-2D regime. 

We stir the cloud with a laser beam which creates 
a repulsive potential with height V^tir ~ &b X 80 nK. 
This is at least four times the local chemical potential 
/ii oc (r) = fi — V(r). The beam has a Gaussian profile 
with a waist of wq — 2.0 (5) ^m, which is larger than the 
local healing length £ = l/*/gn (sa 0.3 yum at the trap 
center), but small compared to the size of the cloud (full 
width at half maximum ps 25 /im) (see Fig. [I]) . We stir 
for typically i s ti r = 0.2 s at constant velocity v in a cir- 
cle of radius r centred on the cloud. The intensity of the 
stirring beam is ramped on and off in « 5 ms without any 



significant additional heating. Once the stirring beam is 
switched off, we let the cloud relax for 0.1 s and measure 
the temperature. 

For each configuration (TV, T, r), we repeat this ex- 
periment for various v from to 2 mm/s. We find two 
different regimes for the response and we show an exam- 
ple of each in Fig.[2j In Fig.[2^i, there is a clear thresh- 
old behaviour with no discernable dissipation below a 
critical velocity. In contrast, in Fig.[2]b, the tempera- 
ture increases without a threshold. We identify these 
behaviours as the superfluid and normal response, re- 
spectively. To model these data we choose for a given 
configuration the fit function 

T(v) = T V=Q + k ■ • max[(« 2 - v 2 c ), 0], (1) 

which describes the heating of a 2D superfluid in the pres- 
ence of a moving point-like defect [H] . In equation 
the three fit parameters are the temperature at zero ve- 
locity T v=0l the heating coefficient k, and the critical 
velocity v c . In the normal state, the fit finds v c ~ 
and the according quadratic heating stems from the lin- 
ear scaling of the drag force. Scattering of photons from 
the stirring beam leads to a 'background heating' of less 
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FIG. 3: Superfluid behaviour across the BKT transi- 
tion, a, The critical velocities v c obtained from the curves 
as in Fig.Q plotted versus the single parameter fiioc/ksT, 
which is the relevant quantity due to the scale invariance of 
the weakly-interacting 2D Bose gas. Our data show an ap- 
parent threshold between critical velocities compatible with 
zero and clearly non-zero critical velocities. It is located at 
Mioc I 'k^T ~ 0.24 (dashed line), somewhat above the predic- 
tion (/iioc/fcBr) c = 0.15 for the BKT transition in an infinite 
uniform system [TS] (the grey shaded area indicates the nor- 
mal state by this prediction). The x error bars indicate the 
region of (J,i oc /ksT that is traced by the stirring beam due 
to its size (using the \j\fe. width of the beam) and due to 
the 'background heating'. The y error bar is the fitting er- 
ror. The inset to a shows the critical velocity plotted versus 
the stirring radius r. Due to the different atom numbers and 
temperatures of the clouds, we can find superfluid or normal 
behaviour for the same radius, b, The heating coefficient n as 
a function of [Hoc/k^T for the normal data (red circles) and 
the superfluid data (blue circles). The red solid line shows a 
fit of k linear in the normal density, as expected from a single- 
particle model. The blue dashed line shows an empirical fit 
quadratic in the superfluid density. The calculation for the 
densities assumes T = 90 nK and the densities are averaged 
over the size of the stirring beam. 



than 10 nK compared to the temperature To before stir- 
ring. In the following, we use the mean temperature 
T = (T v —o + To)/2 to characterize the cloud. 

In Fig.[3j we summarize our data obtained for dif- 
ferent configurations (N, T, r). We show in Fig. [3^. 
the fitted critical velocities versus the single parameter 
fii oc (r) / k^T . The relevance of this parameter results 
from (i) the local character of the excitation, so that the 
response of the fluid to the moving perturbation is ex- 
pected to be similar to that of a uniform gas with the 
same temperature and the chemical potential fj.i ac , (ii) 



the scale invariance of the weakly-interacting 2D Bose 
gas, whose thermodynamical properties do not depend 
separately on /x and T, but only on the ratio [i/k^T 
P31 HS1 HZ]- In particular, this ratio is univocally re- 
lated to the phase space density, and thus characterizes 
the degree of degeneracy of the cloud. 

Quite remarkably, the ensemble of our data when plot- 
ted as a function of ^\o C /k^T shows an apparent thresh- 
old between values compatible with zero and clearly non- 
zero values. This threshold is located at ^loc/k^T ss 0.24, 
somewhat above the prediction (fi/k^T) — 0.15 for the 
superfluid phase transition in a uniform system |15j with 
g = 0.093. The deviation may be attributed to (i) finite- 
size effects in our trapped atomic clouds; (ii) the non- 
zero width of the stirring beam, which extends over a 
significant range of pL\oc/kB,T (horizontal error bars in 
Fig.gi). The emergence of quasi-long range coherence 
in a trapped 2D system was found [6] to agree with the 
prediction for a homogeneous system. If the same is true 
for the superfluid behaviour, this means that effect (ii) 
dominates and that the stirrer has to stand entirely in 
the superfluid core in order to lead to a non-zero critical 
velocity. 

We limit the presented stirring radii to r > 10 /im 
such that the stirring frequencies w — v/r for the rele- 
vant velocities v ~ v c are well below u> r . Indeed, smaller 
radii correspond to a larger centripetal acceleration. This 
could lead to additional heating via the phonon analog of 
synchrotron radiation, as observed in the formally similar 
context of capillary waves generated by a rotating object 

m 

For a homogeneous system, the value of the critical 
velocity is limited by two dissipation mechanisms, the 
excitation of phonons or vortices. For a point-like ob- 
stacle |16j . phonon excitation dominates and v c is equal 
to the speed of sound, given in the T — limit by 
c s = h^/gn/m (rs 1.6mm/s for n = 50 atoms//im 2 ) (this 
situation is described by the celebrated Landau criterion 
[7]). When the obstacle size wo increases and becomes 
comparable to £, dissipation via the nucleation of vortex- 
antivortex pairs (vortex rings in 3D) becomes significant 
[IWI21) . The corresponding v c is then notably reduced 
with respect to c s . In the limit of very large obstacles 
(wq ^> £), an analytical analysis of the superfluid flow 
stability yields v c ~ H/mwo <C c s [53]. With an 
obstacle size wq > £, our experimental situation is in- 
termediate between these two asymptotic regimes. For 
a non-homogeneous system like ours with the stirring 
obstacle close to the border of the expected superfluid 
regime, one can also excite surface modes [2H [2S] , which 
constitute an additional dissipation mechanism. 

Our measured critical velocities are in the range 0.5- 
l.Omm/s, i.e., v c /c s — 0.3 — 0.6. By contrast, previous 
experiments in 3D clouds found lower fractions v c /c s ~ 
0.1 [5]. The difference may be due to the larger size of 
the obstacles that were used, and to the average along the 
axis of the stirring beam of the density distribution in the 
3D gas [55] . The dominant dissipation mechanism could 
be revealed by e.g. directly observing the created vortex 
pairs as in Ref. [12] or interferometrically detecting the 
Ccrenkov-like wave pattern for v > c s as in experiments 
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with a non-equilibrium 2D superfluid of exciton-polariton 
quasi-particles [27]. 

Fig.[3}D shows the fitted heating coefficients k for the 
normal (red circles) and superfluid data (blue circles). 
In the normal region, we expect the heating to scale lin- 
early with the normal density n no |10j . Using the predic- 
tion of [15 for n no (averaged over the size of the stirring 
beam) we fit k = ax ■ h no and obtain a% m 3 • 10 -6 nK • s. 
This value is in reasonable agreement with the predic- 
tion of a model [10] of a single particle with a thermal 
velocity distribution of mean v — v / 7rk B T/2m collid- 
ing with a moving hard wall of width L = wq yielding 
a x = 16mLv/irNk B ~ 6 • 1(T 6 nK • s (for N = 65000 and 
T = 90 nK) . In particular our data nicely reproduce the 
maximum of n no around the expected superfluid transi- 
tion point. In the superfluid case and v > v c , we empir- 
ically fit a quadratic scaling of the heating with density 
k = (22 * ^sf an d nn( i °2 = 8 • 10~ 9 nK • s • /im 2 . In princi- 
ple, one could develop a more refined model to describe 
the superfluid region, by taking into account the coexis- 
tence of the normal and superfluid states via the sum of 
two heating terms. However, within the accuracy of our 



data, we did not find any evidence for the need of such a 
more refined description. 

In conclusion, we have presented a direct proof of the 
superfluid character of a trapped 2D Bose gas. We ob- 
served friction-less motion of an obstacle below a critical 
velocity if the local chemical potential and the temper- 
ature are such that the ratio nioc/ksT is larger than a 
threshold value. An interesting extension of our work 
would be the study of superfluidity from the complemen- 
tary point of view of persistent currents, by adapting to 
2D the pioneering experiments performed in 3D toroidal 
traps [M30]. 
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